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I. INTRODUCTION 

Testing the nucleon structure continues to be one of the 
most challenging tasks for contemporary experiments in 
particle physics. Pion, photon and electron scattering 
off nucleon targets reveal the non-pointlike nature of the 
nucleon by measuring the interactions that take place 
amongst the nucleon's constituents. 

The lowest-lying excited state of the nucleon, the 
A(1232) baryon, plays an equally important role. It is 
produced in such experiments at LEGS, BATES, MAMI 
and Jefferson Lab. Its mass and decay width are experi- 
mentally established. Because the A mainly decays into 
ttN and much less into 7 AT, pionic effects are expected to 
contribute significantly to its properties. The very small 
mean lifetime of the A translates into a highly unstable 
electromagnetic transition A7A, making a measurement 
of its electromagnetic properties very difficult. For exam- 
ple, the Particle Data Group quotes an estimate that " is 
only a rough guess of the range " for the magnetic moment 
Ha++ — 3.7...7.5/Xjv m 7r+ P — > 7r + 7p experiments [1]. 
The MAMI result for the A + magnetic moment, ob- 
tained in pion radiative photoproduction 7^ — > it j'p, 
is given by fi&+ = 2.7^'$ which includes both ex- 
perimental and theoretical errors [2] . Information on A 
and A - static electromagnetic properties is totally miss- 
ing, and there are no experimental results for the evolu- 
tion of the charge and magnetic properties with Q 2 ^ 0. 
On the other hand, further insight has been achieved 
through the measurement of the A^A transition [3-10], 
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where knowledge of the helicity amplitudes, the electric 
quadrupole and the Coulomb quadrupole form factors of 
the transition allows for an extraction of the A's electric 
quadrupole moment [5]. 

Theoretically the description of the A has been very 
challenging as well. Its properties have been studied in 
quark model calculations [11-22], Skyrme models [23- 
25], chiral cloudy bag models [26, 27], as well as chiral 
effective field theory [28-33]. In the absence of (accu- 
rate) experimental information, model predictions can be 
checked by lattice QCD. While the current-quark mass 
dependence of the baryon decuplet's static electromag- 
netic properties and related issues of their chiral ex- 
trapolation have been studied for quite some time [34- 
38], lattice results for the electromagnetic form factors' 
Q 2 — evolution for a wider range of photon momenta have 
become available only recently [39-41] . 

Understanding the structure of the A baryon, its defor- 
mation from sphericity, and the connection to the prop- 
erties of the nucleon via the A7A quadrupole transitions 
must be complemented by extensive research of the elec- 
tromagnetic vertex A7A. Comparative studies of both 
the A/7 A and A7A transition will reveal to which extent 
the deformation of the A-baryon is provided by orbital 
angular-momentum components of its constituents. Nat- 
urally such a study will also shed light onto the nature 
of the A as a pure quark state, rather than a molecular 
state. Finally, perhaps the most important issue to be 
answered is the chiral cloud content of the A-baryon. 

In connection to this, a QCD-motivated quark-core 
analysis of A electromagnetic form factors within the 
framework of Dyson-Schwinger equations (DSEs), to- 
gether with hadronic bound-state equations, is expected 
to provide further insight. Dyson-Schwinger equations 
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constitute a fully self-consistent infinite set of coupled 
integral equations for QCD's Green functions. They 
provide a tool to access both perturbative and non- 
perturbative regimes of QCD; see [42-44] for reviews. 
The most prominent phenomena emerging in the latter 
are dynamical chiral symmetry breaking, confinement, 
and the formation of bound states which require a non- 
perturbative treatment. 

Hadrons and their properties are studied in this ap- 
proach via covariant bound-state equations, see [45-47] 
and references therein. While mesons can be described 
by solutions of the gg-bound-state Bethe-Salpeter equa- 
tion (BSE), the case of a baryon is more involved. The 
three-body equivalent of the BSE is the covariant Fad- 
deev equation. It was recently solved for the nucleon 
mass by implementing a rainbow-ladder (RL) truncation, 
i.e. a dressed gluon-ladder exchange kernel between any 
two quarks, thereby enabling a direct comparison with 
corresponding meson studies [48, 49]. 

In the absence of a solution for the A in this framework, 
a practicable simplification of the problem is based on the 
observation that the attractive nature of quark-antiquark 
correlations in a color-singlet meson is also attractive 
for 3c quark-quark correlations within a color-singlet 
baryon. This provides the tools for studying the three- 
quark problem by means of a covariant quark- 'diquark' 
bound-state BSE [50, 51]. At the current level of com- 
plexity, the importance of meson-cloud effects in the chi- 
ral and low-momentum structure of hadrons is not yet 
accounted for, hence the framework aims at a descrip- 
tion of the hadronic quark core. 

In the present work we adopt this procedure to com- 
pute the electromagnetic properties of the A(1232). This 
augments our previous investigations of quark-core con- 
tributions to the A-baryon mass [52, 53] and nucleon 
mass and form factors [54, 55]. At the same time it rep- 
resents an intermediate step towards the description of 
the nontrivial N — > A7 transition. 

We organize the manuscript as follows: in Section II 
we briefly summarize the Poincare-covariant Faddeev ap- 
proach to baryons and its simplification to a quark- 
diquark picture; and we collect the ingredients of the 
covariant quark-diquark BSE. In Section III we discuss 
the properties and construction of the A electromagnetic 
current operator. In Section IV we present and comment 
on the results for the A electromagnetic form factors and 
static properties. We also compare our results with a 
selection of lattice-QCD results as well as the available 
experimental data for the A and Q baryon. Technical de- 
tails of the calculation are collected in Appendices A-C. 
Throughout this paper we work in Euclidean momentum 
space and use the isospin-symmetric limit m u = rrid- 
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FIG. 1: (Color online) The quark-diquark BSE, Eq. (1) 



II. QUARK-DIQUARK FADDEEV-EQUATION 
FRAMEWORK 

Baryonic bound states correspond to poles in the three- 
quark scattering matrix. The three-quark bound-state 
amplitude is defined as the residue at the pole associ- 
ated to a baryon of mass M. It satisfies a covariant ho- 
mogeneous integral equation, which, upon neglecting ir- 
reducible three-body interactions, leads to the covariant 
Faddeev equation [56] that traces the binding mechanism 
of three quarks in a baryon to its quark-quark correla- 
tions. 

A viable truncation of the Faddeev equation intro- 
duces diquarks as explicit degrees of freedom. It has been 
demonstrated that the same mechanism that binds color- 
singlet mesons is suitable to account for an attraction 
in the corresponding diquark channels within the baryon 
[57, 58]. In particular, a color-singlet baryon emerges as a 
bound state of a color-triplet quark and color-antitriplet 
diquark correlations which are implemented via a sep- 
arable sum of pseudoparticle-pole contributions in the 
quark-quark scattering matrix. 

This procedure leads to a quark-diquark BSE on the 
baryon's mass shell, where the lightest diquarks, i.e. the 
scalar + and axial-vector 1 + ones, have been used to 
describe the nucleon. The spin— 3 /2 and isospin— 3 /2 fla- 
vor symmetric A necessitates only axial-vector diquark 
correlations; its quark-diquark BSE reads [51, 59] 

$""(p, P) = f K^ip, k, P) S(k q ) D»%k d ) $^(fc, P) , 

k 

(1) 

where P is the total baryon momentum, k q , kd are quark 
and diquark momenta, p, k are the quark-diquark rela- 
tive momenta, and J, denotes J <Pk/ (2ir) A . Greek su- 
perscripts represent Lorentz indices, Greek subscripts 
fermion indices. The amplitudes $^o(p, P) arc the 
matrix-valued remainders of the full quark-diquark am- 
plitude $^(p,P)u£(P) for the A, where u£(P) is a 
Rarita-Schwinger spinor describing a free spin-3/2 par- 
ticle with momentum P. 

In order to solve Eq. (1) one needs to specify the 
dressed-quark propagator S, the axial-vector diquark 
propagator D p<y , and the axial- vector diquark ampli- 
tude r" and its charge-conjugate T u which appear in the 
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FIG. 2: (Color online) Quark DSE (4) and diquark BSE (7) in rainbow-ladder truncation. 



quark-diquark kernel: 

K^(p, k, P) = T v {k r , k d ) S T (q) T»(p r ,p d ) , (2) 

where subscripts "r" denote quark-quark relative mo- 
menta and " d" diquark momenta. The mechanism which 
binds the A and is expressed through Eqs. (1-2) is an it- 
erated exchange of roles between the single quark and any 
of the quarks contained in the diquark. This exchange is 
depicted in Fig. 1. 

A solution of the quark-diquark BSE for the A was 
presented in Ref. [52] and corresponding results for its 
mass were reported therein. In the following subsections 
we proceed by recollecting the ingredients of Eq. (1). 

A. Quark propagator and quark-gluon coupling 

The fundamental building block which appears in 
Eqs. (1-2) and connects the quark-diquark model and re- 
sulting hadron properties with the underlying structure 
of QCD is the dressed quark propagator S(p). It is ex- 
pressed in terms of two scalar functions, 

S- 1 (p) = A(p 2 ) (it + M(p 2 )), (3) 

namely the quark wave-function renormalization 1/A(p 2 ) 
and the quark mass function M(p 2 ). Dynamical chiral 
symmetry breaking becomes manifest through a non- 
perturbative enhancement of both dressing functions 
M(p 2 ) and A(p 2 ) at small momenta which indicates the 
dynamical generation of a large constituent-quark mass. 

Such a dynamical enhancement emerges in the solution 
of the quark DSE, cf. Fig. 2: 

s afi(p) = Z 2 (¥ + m ) a p+ J fc aa >p>p(p,q) S a >f3>(q) , (4) 

Q 

where Z 2 is the quark renormalization constant and m 
the bare current-quark mass which constitutes an input 
of the equation. The interaction kernel K. includes the 
dressed gluon propagator as well as one bare and one 
dressed quark-gluon vertex. 

In principle, the dressed gluon propagator and quark- 
gluon vertex could be obtained as solutions of the infinite 
coupled tower of QCD's DSEs, together with all other 



Green functions of the theory. In practical numerical 
studies one employs a truncation: only a subset of the 
infinite system of equations is solved for explicitly; Green 
functions appearing in the subset but not solved for are 
represented by substantiated ansatze. 

In connection with meson properties, e.g. to establish 
the pion as the Goldstone boson of spontaneous chiral 
symmetry breaking, it is imperative to employ a trun- 
cation that preserves the axial-vector Ward-Takahashi 
identity. The latter connects the kernel of the quark DSE 
with that of a meson BSE, ensures a massless pion in the 
chiral limit and leads to a generalized Gell-Mann-Oakes- 
Renner relation [60, 61]. Such a symmetry-preserving 
truncation scheme was described in [62, 63], and its low- 
est order is the rainbow-ladder (RL) truncation which 
amounts to an iterated dressed-gluon exchange between 
quark and antiquark. It has been extensively used in 
Dyson-Schwinger studies of hadrons, see e.g. [64, 65] 
and references therein. The RL truncation retains only 
the vector part ~ 7^ of the dressed quark-gluon vertex. 
Its non-perturbative dressing, together with that of the 
gluon propagator, is absorbed into an effective coupling 
a(k 2 ) which is modeled. The kernel K of both quark DSE 
and meson BSE then reads: 

2 47ra(fc 2 ) v 
KaafW = Z 2 — - 2 — T£ 7 ^ Q , ~fpp, , (5) 

where T£ v = 6^ v — k^k is a transverse projector with 
respect to the gluon momentum k — q — p, and k^ — 
k^ /VP denotes a normalized 4- vector. 

At large gluon momenta, the effective coupling a(k 2 ) 
is constrained by perturbative QCD; in the deep in- 
frared, its behavior is irrelevant for hadronic ground 
states [66]. At small and intermediate momenta it must 
exhibit sufficient strength to allow for dynamical chiral 
symmetry breaking and the dynamical generation of a 
constituent-quark mass scale. We employ the frequently 
used ansatz [67] 

a(k 2 ) = w( 1 pj 2 e- r > 2 (£) +a uv (fc 2 ), (6) 

where the second term reproduces the logarithmic de- 
crease of QCD's perturbative running coupling and van- 
ishes at k 2 = 0. The first term supplies the necessary in- 
frared strength and is characterized by two parameters: 
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FIG. 3: (Color online) Effective coupling a(k 2 ) of Eq. (6), 
evaluated for A = 0.98 GeV corresponding to the u/d— quark 
mass and in the range r) £ [1.6, 2.0]. 



an infrared scale A and a dimensionless width parameter 
r), cf. Fig. 3. (They are related to the infrared parame- 
ters of Ref. [52] via c = (A/A ) 3 and ui = r]~ 1 A/A , with 
A = 1 GeV.) 

The interaction of Eq. (6) provides a reasonable de- 
scription of pseudoscalar-meson, vector-meson and nu- 
cleon ground-state properties if the scale A is adjusted 
to reproduce the experimental pion decay constant and 
kept fixed for all values of the quark mass (see [52, 64, 68] 
and references therein). The corresponding value is A = 
0.72 GeV. Furthermore, these observables have turned 
out to be insensitive to the shape of the coupling in the 
infrared [64, 67]; i.e., to a variation of the parameter r\ 
around the value 77 sa 1.8. 

Other quantities, most notably the masses of axial- 
vector and pseudoscalar isosinglet mesons, are not repro- 
duced so well in a RL truncation. Efforts to go beyond 
RL have been made, and are underway (see e.g. [69-71]), 
but typically require a significant amplification of numer- 
ical effort. In recent studies certain additional structures 
in the quark-gluon vertex, and for consistency also in the 
quark-antiquark kernel, have proven capable to provide 
a better description of such observables as well [70-74] . 

On the other hand, substantial attractive contributions 
come from a pseudoscalar meson cloud which augments 
the 'quark core' of dynamically generated hadron observ- 
ables in the chiral regime, whereas it vanishes with in- 
creasing current-quark mass. A viewpoint explored in 
Ref. [75] was to identify RL with the quark core of chi- 
ral effective field theory which, among other corrections, 
must be subsequently dressed by pion-cloud effects. From 
this perspective a coincidence of RL results in the chiral 
region with experimental or lattice data becomes objec- 
tionable. The properties of a hadronic quark core were 
then mimicked by implementing a current-mass depen- 
dent scale A(m) which is deliberately inflated close to 



the chiral limit, where Awl GeV. As a result, mass- 
dimensionful n, p, N and A observables were shown 
to be consistently overestimated and mostly compatible 
with quark-core estimates from quark models and chi- 
ral perturbation theory; for a detailed discussion, see 
[52, 54, 75]. 

In the present work we employ this 'core model' of 
Ref. [75] to compute the A7A transition properties. 
However, as we will argue in Section IV, the distinction 
between the core model (A(m)) and the fixed-scale ver- 
sion (A = 0.72 GeV) becomes mostly irrelevant once the 
scale of dynamical chiral symmetry breaking (here: the 
mass of the A) has been set and all dimensionful quan- 
tities are expressed in terms of this scale. We finally 
stress that through Eq. (6) all parameters of the interac- 
tion a(k 2 ) are fixed by using information from tt— and 
p— meson core properties only. 



B. Diquarks 

Various theoretical approaches as well as experimen- 
tal observations indicate that the strong attraction be- 
tween two quarks to form diquarks within a baryon is a 
key feature for a better understanding of hadron proper- 
ties [76, 77]. This entails that the quark-quark scatter- 
ing matrix is dominated by diquark degrees of freedom 
at small spacelike and timelike values of the total two- 
quark momentum P. A certain singularity structure in 
the timelike region indicates the presence of diquark mass 
scales within a baryon. In the simplest case such a struc- 
ture can be realized through timelike diquark poles at 
certain values of P 2 , i.e. P 2 — —m 2 c , P 2 = —m 2 v , which 
characterize the lightest diquarks, namely the scalar and 
axial- vector ones. Diquarks carry color and are hence 
not observable; yet such a pole structure does per se not 
contradict diquark confinement, see e.g. [43]. 

In the present context, timelike diquark poles emerge 
as an artifact of the RL truncation which does not per- 
sist beyond RL [63]. In complete analogy to a meson 
or baryon case, it nevertheless allows to derive bound- 
state equations at these poles which determine the on- 
shell scalar and axial-vector diquark amplitudes that ap- 
pear in Eq. (2), together with their masses. The assump- 
tion that this separable structure of the scattering matrix 
persists for all values of P 2 is the underlying condition 
which simplifies the Faddeev equation to a quark-diquark 
model. 

Compared to the nucleon, only an isospin-1 diquark 
can contribute to the isospin- 3 /2 A amplitude which ex- 
cludes the involvement of a scalar diquark. The Bcthc- 
Salpeter equation for the on-shell axial-vector diquark 
amplitude reads 

K (P, P ) = J 1C aa >w{S{q + )T»( q ,P)S T { q _)} aip ,, 
1 

(7) 
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FIG. 4: (Color online) Result for the quark-diquark A am- 
plitude. The plot shows the zeroth Chebyshev moments of 
the dressing functions /j~ in Eq. (9), where the labeling cor- 
responds to Eq. 18 of Ref. [52]. The remaining components 
/3, fs and fs are small and not displayed. 



C. Quark-diquark A amplitudes 

All ingredients of the quark-diquark BSE (1) are now 
specified: the quark propagator as obtained from the 
quark DSE (4), the diquark amplitude as a solution of 
the diquark BSE (7), and the diquark propagator from 
Eq. (8); all obtained within a RL truncation that involves 
the effective coupling a(k 2 ). To compute the amplitude 
and mass of the A baryon numerically, the structure of 
the on-shell quark-diquark amplitude ^ u must be spec- 
ified. It is decomposed into 8 covariant and orthogonal 
basis elements: 

8 

^(p, P)=J2 fk(P 2 ,P ■ P) < p { P , P) P^(P) , (9) 

k—1 

where P is the A onshcll momentum with P 2 = — M A 
and P = P/ (jMa). The Rarita-Schwinger projector onto 
positive energy and spin-3/2 is given by 



P""(ri = A_ ( n(r^-i^^ 



(10) 



where P is the diquark momentum, p is the relative mo- 
mentum between the two quarks in the diquark bound 
state and q± = ±q + P/2 are the quark momenta. The 
equation has the same shape as a vector-meson BSE; in 
a RL truncation the kernel K is given by Eq. (5) . The 
inherent color structure of the kernel leads to prefactors 
4/3 and —2/3 for the integrals in (4) and (7), respec- 
tively. Poincare covariance entails that the axial-vector 
diquark amplitude does not only consist of its dominant 
structure ~ 7 A 'C, but involves 12 momentum-dependent 
basis elements which are self-consistently generated upon 
solving Eq. (7). 

The diquark BSE only specifies the on-shell diquark 



amplitude, i.e. at P 
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Diquarks in a baryon 



are offshell. Information on the off-shell behavior of the 
scattering matrix T can be inferred from its Dyson series, 
schematically written as T = JC + J JCSS T. Reinserting 
the separable ansatz T = T^D^ V T U yields an expression 
for the axial- vector diquark propagator D^ V (P): 

D-l(P) = ml {A 8 M „ +pF flv (P) + Q MV (P)} , (8) 

where and F^v are one- and two-loop integrals in- 
volving diquark amplitudes, the quark propagator and 
the RL kernel K., and A and j3 are related to their on- 
shell values; see [47, 52] for details. Eq. (8) completely 
specifies the diquark propagator from its substructure. 
On the mass shell, it behaves like a transverse particle 



pole: D~i(P 2 



,) = (P 2 + m 2 v ) T£ v , whereas for 



offshell momenta, via Dyson's equation, it picks up non- 
resonant contributions which are implicit in the T-matrix 
but cannot be described by a free-particle propagator. 



where A + (P) = (1 + if*)/ 2 is the positive-energy projec- 
tor, Tp 13 — <5 Q/3 — P a pi 3 is a transverse projector with 
respect to the total momentum, and 7^ = Tp^j 13 are 
transverse 7— matrices. 

The details of the calculation as well as results for the 
A mass were reported in Ref. [52]. A partial- wave anal- 
ysis of the quark-diquark amplitude assigns total quark- 
diquark spin and orbital angular momentum quantum 
numbers to each of the 8 basis elements. Upon per- 
forming the flavor-color traces, a standard procedure to 
solve the quark-diquark BSE involves a Chebyshev ex- 
pansion in the angular variable p ■ P and leads to cou- 
pled one-dimensional eigenvalue equations for the Cheby- 
shev moments of the dressing functions fj^. They match 
the BSE solution at P 2 = — M A , i.e. for an eigenvalue 
A BS e(P 2 = -Ml) = 1. 

The A mass in Ref. [52] was calculated for both model 
versions discussed in Sec. II A. Using a fixed scale A = 
0.72 GeV yields the result Ma = 1.28 GeV which is rea- 
sonably close to the experimental value 1.232 GeV. In 
the core version, where A = 0.98 GeV at the u/d mass, 
the result is Ma = 1.73(5) GeV, where the bracket de- 
notes the sensitivity to the infrared width parameter 77. 
From the perspective of chiral effective field theory, pi- 
onic effects should reduce the A mass by merely ~ 300 
MeV which might indicate the relevance of further di- 
quark channels in describing A properties. On the other 
hand the 'core' A in the present approach is not a reso- 
nance since the A — > iV7r decay channel is not accounted 
for in the quark-diquark kernel, and a corresponding non- 
zero width might impact on its mass as well. 

The result for the A amplitude, as obtained in the 
RL-truncated quark-diquark approach, is dominated by 
an s— wave component in its rest frame, cf. Fig. 4. The 
subleading s— , p— and d— wave amplitude components 
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are significantly suppressed compared to this structure 
which corresponds to rf p = S pp in Eq. (9) . A similar 



observation holds for the nucleon amplitude and might „n c n c n 

. ,. , , , ., . . * , ,. . ^dim _ e ^E 2 ^dim _ e u Mi ridim _ e U M 3 

indicate that orbital angular-momentum correlations m { ^ r E 2 ~ ^2 > tT Afi — ~2~M — ' Ms ~ 2M 3 

these baryons' amplitudes are dominated by pionic effects v v x 
which are absent in our setup. 
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III. DELTA ELECTROMAGNETIC FORM 
FACTORS 

A. Electromagnetic current operator 

Having numerically calculated the A-baryon ampli- 
tudes, we proceed with the construction of the A elec- 
tromagnetic current. It can be written in the form 



J p ' pr7 (P,Q) = iF pa (P f ) 



2Ma 



4M A 



F, 



2M A 



§a 
(11) 



which is derived in App. B2. The exchanged photon 
momentum is denoted by Q — Pf — Pi, where Pj and 
Pf are the initial and final momenta of the A and P = 
(Pi + Pf)/2 is its average total momentum. The Rarita- 
Schwinger projectors were defined in Eq. (10). 

The electromagnetic current is expressed in terms of 
four form factors F*(Q 2 ). The experimentally mea- 
sured A form factors - Coulomb monopole Geo, mag- 
netic dipole Gmi, electric quadrupole Ge2, and magnetic 
octupole Gm3 - can be expressed through linear combi- 
nations of the F*(Q 2 ) [78, 79]: 



1 + y) " tF$) - T -{1 + t) (F* - tFZ) 



G 



Mi 



1 + y) (f? + Ft) - f (1 +r) (F* + Ft) 
1 



Ge 2 := {Ft - rF*) - - (1 + r) (F 3 * - rF 4 *) , 
G M3 := (*? + F*) - \ (1 + r) (F 3 * + F 4 *) ■ (12) 
Their static dimensionless values are given by 
G Eo (0) = e A , G Ml (0) = Ma > 



G Ba (0) = Q, G M3 (0) = O. 



(13) 



where e A € {2, 1, 0, —1} is the A charge, /i A its magnetic 
dipole moment, Q the electric quadrupole moment, and 
O the magnetic octupole moment. Equivalently, one has 



F£{0) = e A , 



F${0) = /i A e A 



F 3 *(0)=e A -g, F 4 *(0) =^ A -e A + Q-G. 



(14) 



These form factors are dimensionless. Their dimensionful 
values are given by 



B. Construction of the electromagnetic current 

To compute the electromagnetic properties of the A- 
baryon in a given framework, one must specify how the 
photon couples to its constituents. In the quark-diquark 
context this amounts to resolving the coupling of the pho- 
ton to the dressed quark, to the diquark, and to the inter- 
action between them, where the incoming and outgoing 
baryon states are described by the quark-diquark ampli- 
tudes of Eq. (9). 

The construction of this current is based on a proce- 
dure which automatically satisfies electromagnetic gauge 
invariance [80, 81]. The corresponding diagrams are de- 
picted in Fig. 5 and worked out in detail in App. C. The 
upper left diagram describes the impulse-approximation 
coupling of the photon to the dressed quark and involves 
the quark-photon vertex. The lower left diagram is the 
respective coupling to the diquark and depends on the 
axial-vector diquark-photon vertex. The upper right di- 
agram depicts the photon's coupling to the exchanged 
quark in the quark-diquark kernel, and the lower two 
diagrams its coupling to the diquark amplitudes which 
involve seagull vertices. 

At the level of the constituents, electromagnetic cur- 
rent conservation Q p -J p ' f " T = translates to Ward- 
Takahashi identities which constrain these vertices and 
relate them to the previously determined quark and di- 
quark propagators and diquark amplitudes. Neverthe- 
less, the vertices may involve parts transverse to the 
photon momentum which are not constrained by current 
conservation and yet encode important physics. A self- 
consistent determination of such transverse parts is in 
principle possible but requires certain numerical effort. 
For instance, the quark-photon vertex can be computed 
from its rainbow-ladder truncated inhomogeneous Bethe- 
Salpeter equation which unambiguously fixes its trans- 
verse contribution [82]. As expected from vector- meson 
dominance models, the latter exhibits a p— meson pole at 
Q 2 = -m 2 p . 

In the present calculation we construct the quark- 
photon vertex from its component fixed by the WTI, 
i.e. the Ball-Chiu vertex, augmented by a transverse p- 
meson pole contribution that is modeled after the result 
in [82]. An analogous construction is used for the axial- 
vector seagull vertex. Having fixed those, the axial- vector 
diquark-photon vertex is completely specified. The de- 
tails of the construction are presented in Apps. C 2-C 4. 
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FIG. 5: (Color online) A baryon's electromagnetic current in the quark-diquark model. 
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0.14 
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0.77(3) 
0.88(5) 


2.78(12) 
2.12(8) 


0.70(6) 
0.38(2) 



TABLE I: Results for the A mass and static electromagnetic properties at two different pion masses, where A is the input scale 
of Eq. (6) and the 77 dependence is indicated in parentheses, m^, A and Ma are given in GeV, /za in nuclear magnetons and 
2 The Gi(0) are dimensionless and the Ri are given in (GeVfm) 2 . 
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IV. RESULTS AND DISCUSSION 

We have calculated the electromagnetic properties of 
the A within the decomposition of Fig. 5 for the electro- 
magnetic current. This is achieved by specifying expres- 
sions for the quark-photon vertex, the diquark-photon 
vertices and the seagull terms. All necessary ingredi- 
ents that enter this calculation, such as the quark and 
diquark propagators and the axial-vector diquark and 
A amplitudes, were explained in Section II. Through 
the decomposition of Fig. 5 all resulting A properties 
are traced back to the effective quark-gluon coupling in 
Eq. (6). Specifically, one can investigate the impact of the 
infrared properties, controlled by the width parameter 77 
(cf. Fig. 3), on resulting observables. This is indicated 
by the colored bands in Figs. 6 and 7. 

In Fig. 6 we depict the core contributions to the A + 
electromagnetic form factors, calculated at the physical 
point TOtt = 140 MeV and compared to lattice data at 



three different pion masses [40]. As can be seen from 
the figure, we find a good overall agreement between our 
results and those obtained on the lattice. Because we as- 
sume isospin symmetry, the A ++ , A and A~ form fac- 
tors are simply obtained by multiplying those of the A + 
with the appropriate charges, cf. App. C 5. We note that 
the A in our core calculation is a stable bound state and 
not a resonance, i.e. it does not develop a width. Non- 
analyticities associated with the decay channel A — > Nir 
which would appear for Ma — Mjv > , corresponding 
to the domain < 300 MeV, are therefore absent. The 
same is true for the available lattice results which are 
obtained at pion masses above this threshold. 

Since the A is highly unstable, the available exper- 
imental data for its electromagnetic decays are rather 
poor. Only the magnetic moments of A + and A ++ 
are experimentally known, albeit with large errors. The 
corresponding values for Gmi(0) are 7.3 ± 2.5 (A ++ ) 
and 3.5±£j (A+). Our result G M1 (0) = 3.64(16) com- 
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FIG. 6: (Color online) Electromagnetic form factors of the A. The bands represent the sensitivity to a variation of r\ = 1.8±0.2. 
The results are compared to unquenched lattice data of Ref. [40] at three different pion masses. 



pares well with quark-model predictions and chirally ex- 
trapolated lattice results which typically quote values 
Gmi(O) ~ 3 ... 4, see [22] and references therein. 

The deformation of the A is encoded in its electric 
quadrupole and magnetic octupole form factors. Non- 
relativistically, a negative sign for the electric quadrupole 
moment Gei^S*) indicates an oblate charge distribution 
for the A in the Breit frame. (A different interpreta- 
tion arises in the infinite- momentum frame, see [40]). 
From the measurement of the NjA transition, one can 
infer the value G B2 (0) = -1.87(8) in the large-iV c 
limit [40, 83]; comparable values are predicted by a range 
of constituent-quark models [20]. The impact of pionic 
corrections upon the low-Q 2 behavior of Ge2 is however 
unclear. While the lattice data of Ref. [40] are limited 
by large statistical errors which prevent an accurate ex- 
trapolation to Q 2 = 0, they indicate a negative value for 
G_E2(0) as well. Our result for the electric quadrupole 
moment, Ge2(0) = —1.32(16), is negative and compati- 
ble with these results. We note that Ge2(Q 2 ) develops a 
zero-crossing at Q 2 jM\ ~ 0.6, a feature which is unex- 



pected but not clearly excluded from the available lattice 
results. 

The lattice signal for the magnetic octupole form fac- 
tor Gm3(Q 2 ) is weak and plagued by large error bars, 
especially at low Q 2 . At next-to- leading order in a 
chiral expansion, the magnetic octupole moment van- 
ishes [84]. Our calculation yields a small and negative 
form factor; the corresponding magnetic octupole mo- 
ment is G M3 (0) = -0.26(4). 

An ambiguity arises when comparing our form factor 
results, calculated with an implicit 'core' mass Ma > 
M^ xp , to experimental and lattice data. Form factors 
are dimensionless, hence they can only depend on di- 
mensionless variables. To account for this, the electro- 
magnetic form factors in Fig. 6 are plotted as a func- 
tion of the dimensionless variable Q 2 /M A for our data 
and Q 2 / (M 1 ^) 2 for the lattice data. From another point 
of view, Ma defines an effective scale of dynamical chi- 
ral symmetry breaking. Once such a scale is set (e.g. 
by having numerically computed Ma), all dimensionful 
quantities can be related to this scale. On the lattice, a 
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FIG. 7: (Color online) Static electromagnetic properties of the A. Upper left panel: dimensionless form factors at Q 2 = 0. 
Upper right panel: dimensionless squared radii Ri, given in (GeVfm) 2 . Lower left panel: A magnetic moment, expressed in 
static nuclear magnetons and compared to lattice data [36-38, 40]. Lower right panel: squared electric charge radius compared 
to lattice results [38, 40]. Stars denote experimental values [1]: for the A we plot ^t A ++/2; for the Q~ we show |/Uq-[. The 
dashed vertical lines indicate the positions of the u/d— and strange-quark mass. Note that there is no ss pseudoscalar meson 
in nature; the value m ss = 0.69 GeV corresponds to a meson-BSE solution at a strange-quark mass m s = 150 MeV [61]. 



scale must be denned as well to convert dimensionless lat- 
tice results into physical units. Hence, an unambiguous 
comparison between different theoretical approaches and 
experiment should ideally involve dimensionless quanti- 
ties. 



creasingly dominated by its core. A similar behavior can 
also be observed for Ge2{0) and G^(0) which are neg- 
ative throughout the current-mass range. 

The A magnetic moment is given by 



In the upper left panel of Fig. 7 we show our results 
for the dimensionless form factors Gmi, Ge2 and Gez at 
vanishing photon momentum transfer Q 2 — 0. For two 
quark flavors with equal masses the A - becomes identi- 
cal to the 17~-baryon when evaluated at the value for the 
strange-quark mass, as the Q,~ is a pure sss state. Ex- 
perimentally, Mq = 1.672 GeV and fiQ- = — 2.02(5)/!^ 
which implies |G M i(0)| = M n /M N |/x n - | = 3.59(9). Our 
result for Gmi (0) is almost independent of the current- 
quark mass and agrees reasonably well with the experi- 
mental value for the Q~ . This is an encouraging result 
as pionic effects should have diminished in the vicinity 
of the strange-quark mass above which the baryon is in- 



dim 
Ma 



2M A 



GmM 



2M C * P 



G 



Ml 



(0)^- 



(15) 



where Ma is running with the current-quark mass; hence 
its value in static nuclear magnetons is given by the 
bracket in Eq. (15). To compare with experiment or lat- 
tice, we must again bear in mind that our calculated A 
mass is different from the one which is measured (or ob- 
tained in lattice calculations) , and that the unambiguous 
comparison of magnetic moments is that of the dimen- 
sionless value Gmi(0). To account for this, we plot //a 
in Fig. 7 by replacing Ma in Eq. (15) by the following 
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reference mass: 

fRef / 2 \2 



3m w 



(l + /(mj)), (16) 



with /(m£) = 0.77/(1 + (m^/O^Gey) 4 ) and M = 1.2 
GeV. Eq. (16) reproduces the experimental A and 
masses at = 0.14 GeV and — 0.69 GeV, re- 
spectively, and provides a reasonable representation of 
the dynamical lattice results of Ref. [40] for At the 
u/d— quark mass, our calculated value for the magnetic 
moment is thus /j,a = 2.78(12)/zjy. Because of the con- 
stancy of Gmi(0) it decreases with an inverse power of 
the mass M^ cf (m^). 

The charge radius corresponding to a form factor 
Gi(Q 2 ) is defined as 



6 dd 
'Gi(O) dQ 2 



(he) 



(17) 



2 =o 



where he = 0.197 GeV fm. We display the dimension- 



less squared charge radii Ri 



' in the upper right 



panel of Fig. 7. Again their current-mass dependence is 
found to be weak. The dependence on the infrared width 
r\ arises from which has a sizeable dependence on rj, 
cf. Table I. In the lower right panel of Fig. 7 we display 
the rescaled squared charge radius 



Mi 



Reo 



(18) 



It shows a satisfactory agreement with the lattice data 
at larger pion masses, where the A is mainly described 
by its core. At the u/d mass, its value is 0.70(6) fm 2 . 

It is noteworthy that an investigation of dimensionless 
properties effectively removes the model dependence be- 
tween the 'core' version and the fixed-scale version of the 
coupling in Eq. (6) which was discussed in Section II A. 
In the chiral limit, in the absence of a current-quark or 
pion mass, the input scale A of Eq. (6) is the only relevant 
scale in all our calculations. For instance, hadron masses 
scale with A. Since the properties of a quark core were 
modeled solely by increasing A, the discrepancy between 
the two models disappears when comparing dimension- 
less values. Nevertheless, since both models describe a 
quark core (pionic effects are absent in either of them), 
characteristic chiral features induced by pionic correc- 
tions, such as a logarithmically divergent A charge ra- 
dius, can naturally not be reproduced. 

On a related note, the estimation of chiral cloud ef- 
fects from a pure quark-core calculation is not unam- 
biguous. For instance, adding pionic effects in the quark- 
diquark kernel would affect the Bethe-Salpeter normal- 
ization of the A amplitude, a condition which is equiva- 
lent to charge conservation: Gso(0) = 1. Charge conser- 
vation must always be satisfied, irrespective of whether 
pionic effects are included or not. Hence the effective 
core contribution to the form factors of a A-baryon that 
is dressed by pionic degrees of freedom would be smaller. 



The resulting 'renormalization' of all form factors would 
be independent of Q 2 but would depend on the current- 
quark mass: at large masses, the difference should dimin- 
ish as pionic effects become small. It is reassuring that 
the A magnetic moment and electric charge radius dis- 
played in Fig. 7 agree quite well with experimental and 
lattice results in this domain. 

We have established that the static limits of the four 
form factors depend only weakly on the current-quark 
mass. It is remarkable that the same can also be said for 
their Q 2 — dependence: when plotted over Q 2 JM\, the 
overall shape of Fig. 6 persists throughout the current- 
quark mass range. For instance, the position of the zero 
crossing in Ge2(Q 2 ), expressed through Q 2 /M%, is al- 
most independent of the value of m n . Thus, our re- 
sulting form factors for the Q~ are not materially dif- 
ferent from those of the A once their inherent mass de- 
pendence is scaled out. The same pattern appears to 
be valid for the lattice results as well: when plotted over 
Q 2 /M^, the results at different pion masses follow a sim- 
ilar Q 2 — evolution. 

We finally note that the current setup prohibits an in- 
vestigation of form factors beyond Q 2 /M\ > 2 because of 
timelike and/or complex singularities that appear in the 
quark propagator and quark-quark scattering matrix. Ir- 
respective of its details which are truncation-dependent, 
such a singularity structure is an inevitable feature of 
these Green functions and must be accounted for in order 
to investigate the large-Q 2 domain. As an intermediate 
step it might be worthwhile to study this regime with the 
help of pole-free model propagators, a procedure which 
has been adopted in Ref. [85] in the context of nucleon 
electromagnetic form factors. 



V. CONCLUSIONS 

We provided a calculation of A and the fi electro- 
magnetic form factors in a Poincare-covariant quark- 
diquark approach. All quark and diquark ingredients 
were determined self-consistently from Dyson-Schwinger 
and Bethe-Salpeter equations and thereby related to the 
fundamental quantities in QCD. We employed a rainbow- 
ladder truncation which corresponds to a dressed gluon 
exchange between the quarks inside the diquark. Since 
pion-cloud effects are not implemented, our results de- 
scribe the A quark core. 

The Q 2 — evolution of the electromagnetic form factors 
agrees well with results from lattice QCD. In particu- 
lar, the electric quadrupole and magnetic octupole form 
factors are negative throughout the current-quark mass 
range which, in the traditional interpretation, indicates 
an oblate deformation of the A's charge and magnetiza- 
tion distributions. The dimensionless form factors at van- 
ishing photon-momentum transfer are almost indepen- 
dent of the current-quark mass. At larger quark masses, 
where pionic effects do not contribute anymore, our re- 
sults for the A magnetic moment and electric charge ra- 
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dius are close to those obtained in lattice QCD. More- 
over, the magnetic moment evaluated at the strange- 
quark mass agrees well with the experimental value for 
the 

In summary, the results collected herein show good 
agreement with lattice observations, quark model anal- 
yses and the available experimental data. Near-future 
measurements at MAMI and JLab facilities remain to 
validate these predictions for the A's electromagnetic 
properties. From this perspective, while the N — > A7 
transition is accurately measured, a forthcoming exten- 
sion of our approach to the investigation of spin- 3 /2 — > 
spin-Y2 electromagnetic transitions will constitute a non- 
trivial test for the Poincare-covariant framework pre- 
sented herein. Moreover, a study of the evolution of 
the A electromagnetic form factors at large values of the 
photon momentum transfer is desirable. Naturally this 
constitutes a very challenging task, in particular since 
data collection for electromagnetic properties of baryon 
resonances at large Q 2 is due to begin. 

Finally, in view of a better understanding of the struc- 
ture of baryons, our approach can be methodically im- 
proved by eliminating the diquark ansatz in support of a 
fully Poincare-covariant solution of the three-quark Fad- 
deev equation. We have recently reported on numerical 
results for the nucleon mass in such a framework [48] . In 
parallel we intend to augment our analysis by develop- 
ing compatible tools to incorporate missing chiral cloud 
effects. 
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where a k are the three Pauli matrices. The charge con- 
jugation matrix is given by 



4„,2 



C 1 



c 1 = c~ 



-c, 
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and the charge conjugates for (pseudo-) scalar, (axial-) 
vector and tensor amplitudes are defined as 

T(p,P) :=CT(-p,-P) T C T , 
f a (p,P) :=-CT a (-p,-P) T C T , (A4) 
T 0a (p,P) :=CT a P(~p,-P) T C T , 

where T denotes a Dirac transpose. Four-momenta are 
conveniently expressed through hyperspherical coordi- 
nates: 
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Appendix A: Euclidean conventions 

We work in Euclidean momentum space with the fol- 
lowing conventions: 

4 

p- q = ^2pkqk, p 2 =p-p, i> = p- n j. (Ai) 

A vector p is spacelike if p 2 > and timelike if p 2 < 0. 
The hermitian 7— matrices 7^ = (7^)^ satisfy the anti- 
commutation relations {7^,7^} = 2(S A " y , and we define 

CT ^ = _* [7 M 7l) 7 s = _ 7 V 7 V ■ (A2) 



Appendix B: A Electromagnetic current 

1. General properties 

The matrix-valued electromagnetic current of the A 
can be written in the most general form as 

J^' pCT (P,g) ^T pa {P f )T^ a!3 {P,Q)¥ 3a {P l ). (Bl) 

The vertex involves two momenta, expressed through the 
ingoing and outgoing momenta Pi, Pf or by the average 
momentum P = (Pi + Pf)/2 and photon momentum Q = 
P f - Pi. Since the particle is on-shell, Pf = Pf = -Af|, 
one has 

P 2 = -M 2 a {1 + t), P-Q = 0, (B2) 

where r := Q 2 /(4M^); hence the Lorentz-invariant 
form factors which constitute the vertex can only de- 
pend on the photon momentum-transfer Q 2 . The Rarita- 
Schwingcr projector W al3 (K) for a general momentum K 
is defined via 

P^(J0 := A+(K) (rf - I 7t7^ , (B3) 
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with the transverse projector T^P = <5 Q/3 — K a K@, 
the transverse gamma matrices 7 ° = T^ /3 7 ' 3 , and the 
positive-energy projector A + (if) = (1 + K de- 

notes a normalized 4- vector. The projectors A + (K) and 
P a/3 (JT) satisfy the relations 

A + (K) ft = $ A + (K) = A+(K) , 

P Q/3 (K) K p = K a V afi (K) = , (B4) 

P«0(jjQ 7 /» =7 a P a ^(Jf) =0. 

By contracting Eq. (Bl) with Rarita-Schwinger spinors 
u p (Pf, Sf), u a (Pi, Si) which are eigenstates of P Q/9 : 

V af, (K)u f, (K,s)=u a (K,s), s e {±3/ 2 , ±1/2}, (B5) 

we obtain the current-matrix element (Pf, Sf | J M | Pi, St). 

The positive-energy projector is invariant under charge 
conjugation: 

A + (K):=CA + (-K) T C T = A + (K), (B6) 

since C C T = 1 and C = — 7 M C. The same is true 
for the charge-conjugated Rarita-Schwinger projector: 



and 



P^(AT): = CT a P{-K) T C T = 



T^- 1 3 4^)A + (K) 



= A + (K)[T^~-4^)=V^(K). 



(B7) 



A + (K) commutes with the bracket since 7^ A + (K) = 
A_(Jf)7£ and hence ^^A+{K) = A + (K)^^, 

where A-(K) = (1 — $T)/2 is the negative-energy pro- 
jector. These relations will be useful in the following 
subsection. 



2. Derivation of the electromagnetic current 

The current J^ CT (P,Q) in Eq. (Bl) consists of four 
form factors. T ,i a ^(P,Q) is a five-point function: it has 
three vector legs and two spinor legs, and it depends on 
two momenta. Such a vertex has 144 tensor structures 
of definite parity which are given by all combinations of 
the 36 elements 
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(B8) 



Q^P a QP 6i* a Q 13 

QPQapP 6^Q a 

pppapf) QPpapfi Stl ap0 [ > 

ppQapP §P0pa 

pppaqji §a^pp 



with the four basic spinor structures 
{1, f,Q, [f,<Q]}. 



(BIO) 



Applying the Rarita-Schwinger projectors reduces this 
set to 8 basis elements. First, only the unit matrix 1 sur- 
vives upon sandwiching Eq. (BIO) between the positive- 
energy projectors, since for instance 

A f + f A\ = A f + A + = iMA + A + ' 

A{<$A i + =A{ (%-^ i )Ai = 0, 

where we have abbreviated A l _r = A + (Pij). Moreover, a 
Rarita-Schwinger contraction of those basis elements in 
Eq. (B8) which involve instances of 7 Q or 7^ yields zero 
due to Eq. (B4), while a contraction with Q a or Q@ can 
be reduced to a contraction with P a or P@ since 

(p 1 p \a pa 
papa^p.) = I tf r S) p»g(p.) ^ _Jlp«0(p.) = 

Pa n a \ 

^_ + y_ j p^ ( p) 

^ p«p^(p.) = 91 p^(p) . 

Discarding all basis elements involving 7 Q , Q a or QP 
one then arrives at the following 8 tensor structures: 



ppfia/3 pHQaQP 

Q>*6 a P Q fl Q a Q 13 



6^Q a . 



(Bll) 



A further reduction can be achieved by imposing 
charge-conjugation invariance of the current J A1 ' pcr (P, Q) 
which, because the Rarita-Schwinger projectors are 
charge-conjugation invariant, translates into charge- 
conjugation invariance of the vertex r M,Q ^: 

r>^ a (P,Q) := -Cr"'^(-P,-Q) T C T = T^ a {P,-Q) . 
It reduces the 8 components to 5: 

p^ap P » Q a Q P S^QP-SrtQ". (B12) 

The resulting current is now automatically conserved, i.e. 
qpjp,p<j _ q r p ae fjf^jj tensor structure can be related 
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to the first three by the following identity [78]: 
■A{ [S^QP -S^Q a ] A\ = 



i 

2M 



Ai 



A 



+ • 

(B13) 



As a consequence, the general electromagnetic current 
of the A depends on four form factors F*(Q 2 ) and can 
be written as: 

J»' pa (P, Q) = P pa (P/) r^' Q/5 (P, Q) P /3<T (P J ) , (B14) 

with 

T^(P, Q) = ((F* + F*) - F* V* 

A P p \ Q a Q p 



Using the Gordon identity 



A 1 , = 



(B15) 



(B16) 



M 2M 

finally leads to the expression given in Eq. (11). 

3. Extraction of the form factors 

The electromagnetic form factors are extracted from 
Lorentz- invariant traces of the current JP> a P, Each in- 
dex /x, a, (3 can be contracted with the momenta P, Q 
or a gamma matrix. For the Rarita-Schwinger indices 
a, /3 it is sufficient to consider the contraction with the 
momentum P due to the relations 

P Q P Q ' 9 (P) = ^ P Q/3 (P) , 7 Q P Q/3 (fc) = . (B17) 

The vector index /i can be contracted with either P M and 
7^ ; a contraction with Q' 1 yields zero because of current 
conservation. In addition, all indices can be contracted 
among themselves. Upon performing the Dirac traces 
one arrives at the following Lorentz-invariant scalars: 



si 

S3 
Si 



= Tr { J^' al3 } pvp a pi 3 
= Tr {J^ aa }P^ 

= Tr {J^ afs ^} P a P f) 
= Tr{J^ QQ 7 £} . 



(B18) 



Further possibilities such as Tr{ P }P^p a pP or 
Tr {J^ Q }P Q = Tr{J^' Q ^}P Q linearly depend on 
(B18). We used P = P/(iM) to arrive at dimensionless 



quantities and 7^ := Tp"j" to project out the compo- 
nent in P— direction. A computation in the Breit frame, 
where 



P = iM N /l + re 4 
Q = 2Aiy?e 3 



P = e 4 . 
Q = e 3 . 



(B19) 



yields the result: 
4zr 
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Since the S; are Lorentz-invariant, this result is frame- 
independent. Its inversion yields for the electric form 
factors 



s 2 - 2si 



Gr, — , 

4zyr+7 
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8iT 2 yTT7 
and for the magnetic form factors: 
9i 

Gmi = 407 (S4 ~ 2s3) ' 
Gm3 " 16^ 



TS 2 



(B20) 
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Appendix C: Diagrams in the quark-diquark model 
1. General decomposition 

In this appendix we collect the ingredients of the A 
electromagnetic current operator in the quark-diquark 
model which are depicted in Fig. 5. The explicit form of 
the current is given by a sum of impulse-approximation 
diagrams (left panel in Fig. 5) and two-loop contribu- 
tions which represent the photon's coupling to the quark- 
diquark kernel (right panel of Fig. 5): 

j^ pa = J $ pa (Pf,Pf) (x q + x dq r^ ^{pi, p^ + 

$ pa (p f ,P f )X£ af} ^( Pl ,p). (Cl) 



Here, Pi and Pf = Pi + Q are incoming and outgo- 
ing on-shell A momenta. The relative momenta pi and 
Pf are independent loop momenta in the two-loop dia- 
grams; in the one-loop diagrams they are related to each 
other: Pf — Pi = (1 — £) Q for the quark diagram and 
Pf —Pi = Q for the diquark diagram, where £ £ [0, 1] 
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is an arbitrary momentum-partitioning parameter which 
must be specified prior to solving the BSE. a,j3= 1 ... 4 
are the diquark's Lorentz indices. The quark-diquark 
amplitudes $ pQ are the solutions of the quark-diquark 
Bethe-Salpeter equation (1). The ingredients of Eq. (CI) 
are given by 



= S(p + )TZ(p + ,p_)S(p^D^(k_) 



(C2) 



X 



dq 



X£ ap = D aa '(k + ) S(p+) K^ a ' f3 ' S(p^)D^ l3 {k_) . 

and depend on the quark-photon vertex and the axial- 
vector diquark photon vertex r^ a/3 . The quark and di- 
quark momenta are: 



P- = Pi +£Pi, 
P+ =Pf+£ p f> 



k- = -Pi + (1 - Pi , 
k+ = -p f + (l-£)P f 



The gauged kernel K^ ,a ^ contains the exchange-quark 
diagram and the seagull vertex Af' a : 



K^ afi = (K EX + K SG + K 



SG' 



(C3) 



with the quark propagator's dressing functions A(p 2 ) and 
B{p 2 ) = M(p 2 )A(p 2 ). The first part is the Ball-Chiu 
vertex [86]. The transverse contribution T p is not con- 
strained by the WTI except for the fact that it must van- 
ish at Q 2 = 0. In the present work we use a phenomeno- 
logical ansatz which is modeled after a RL-truncated 
inhomogeneous BSE solution for the quark-photon ver- 
tex [82] . The latter includes a self-consistently generated 
vector-meson pole at Q 2 — —m 2 whose contribution sig- 
nificantly increases the charge radii of pseudoscalar and 
vector mesons [82, 87]. The ansatz reads 



1 T 



9 P x 



1 



e -9(x) rpllU 



(C7) 



with 



where r^ c is the p— meson amplitude as obtained from its 
homogeneous Bethe-Salpeter equation, x = Q 2 /m 2 , and 

g p = v2 rrip/fp with the computed p mass and decay con- 
stant. An additional function e~ 9 ^ was implemented to 
optimize agreement with the vertex-BSE solution at low 
Q 2 and nucleon form factor phenomenology at interme- 
diate Q 2 , see [54] for details. 



3. Seagulls 



= I*(r + , [S(q+) r> + , q-)S{q-)\ r Q (r_, k+) 

K$g p = (r' + ,k-,Q) S( q+ ) T T<* (r_ ,k+), 

K w0 = T^(r + , fc_) S(q-) T M^' a (r'_, k + , Q) , (C4) 

and momenta: 

p±-q T , _P±~Q± 



Q± 



Pt 



For explicit calculations we work in the Breit frame (B19) 
where ingoing and outgoing nucleon have opposite 3- 
momenta and the photon consequentially carries zero en- 
ergy. 



The seagull vertices M M ' Q satisfy a Ward-Takahashi 
identity similar to Eq. (C5) which involves differences of 
axial- vector diquark amplitudes [65, 81, 88]. The result- 
ing vertex as derived from this WTI has again the form 



(C8) 



In a similar spirit as before we include a phenomeno- 
logical transverse p— meson part which optimizes 

agreement with nucleon form factors at intermediate Q 2 
but is irrelevant on the domain Q 2 < 1 GeV 2 . The de- 
tailed structure and derivation of (C8) can be found in 
App. A.8 of Ref. [47]. 



2. Quark-photon vertex 

The general expression for the quark-photon vertex 
(fc + , fc_) is derived from the Ward-Takahashi identity 

QT£(fc+,fc_) = S-\k + )- S-\k^) (C5) 

and by imposing regularity at Q 2 = 0. It is given by 

r£(fc+,A_)=*7^E A + 2^(#A A + A B )+f£ 1 (C6) 

where k + and fc_ are outgoing and incoming quark mo- 
menta, k = (fc + +fc_)/2, Q = k + — fe_, and 

_ F(k 2 + ) + F(k 2 _) _ F{k\)-F{k 2 _) 
Ef: " 2 ' Af: " kl-kl ' 



4. Diquark-photon vertex 

The procedure which allows to derive the electromag- 
netic current for a composite particle can be generalized 
for an arbitrary composite vertex at off-shell momenta. 
Once the propagator of the composite object is known, 
the vertex is constructed as the 'gauged' inverse propa- 
gator, i.e. such that the photon couples to all its con- 
stituents. The Ward-Takahashi identity for the vertex is 
then automatically satisfied [80, 89]. 

Starting from the expression for the axial-vector di- 
quark propagator, Eq. (8), the axial- vector diquark- 
photon vertex is constructed by coupling the photon to 
each of the constituents of the one- and two-loop integrals 
F^y and Q^. These involve the impulse-approximation 
coupling to the quark, specified by the quark-photon ver- 
tex, and the coupling to the diquark amplitudes, i.e. 
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the seagulls. The seagull terms vanish on the diquark's 
mass shell — hence the RL-consistent meson- or diquark- 
photon current is described by the impulse approxima- 
tion — but must be considered for off-shell momenta in 
order to satisfy the WTI. 

Having determined the quark-photon and seagull ver- 
tices previously, the diquark-photon vertex is completely 
specified. It is given in App. A. 7 of Ref. [47]. 



5. Color, flavor and charge coefficients 

The current matrix diagrams of App. C 1 still have 
to be equipped with color and flavor-charge coefficients. 
The full Dirac, color and flavor structure of the A quark- 
diquark amplitude is given by 



Sab 
V3 



L ab I 



(C9) 



where is the Dirac-momentum part of Eq. (9) and 
Sab its color factor, with A,B = 1,2,3 corresponding 
to the quark and diquark legs. The flavor matrices t!, b , 
where the index i = 1,2,3 belongs to the axial- vector 
diquark's three symmetric isospin-1 states, the index a 
1 . 2 to the quark's two isospin- 1 ^ states and the index 
b = 1 ... 4 to the four isospin- 3 /2 states of the A, are 
given by 




i 

\/3 








t z := 




73 
1 



Denoting the four isospin- 3 /2 unit vectors e r by 



e 1 1 = 



e T = 



e = 




e = 




the isospin-Ya basis by 



u = 



d = 



and the isospin-triplet diquark matrices by 



s 1 = uu^ 



s 2 = -= (iid* + du 1 ) , s 3 = dd t , 
v2 



the contraction t 1 e r represents the remainders upon re- 
moving the diquark flavor matrices s 1 from the full three- 
quark flavor wave function Xa=i s 1 ® t 1 e r . The corre- 
sponding 'super'-vectors A ++ , A + , A , A~ whose entries 
are (A r )' = t 1 e r are then given by 







( o ^ 






\/t u 


i 


\/f d 




! : ) 


V J 




K 73 u J 




w 



These definitions allow to write down concise expres- 
sions for the color and flavor traces of the electromagnetic 
current. The color traces for the impulse approximation 
and exchange/seagull diagrams are given by 



Sba 



Sab 



1 



Sba £aed £ceb Scd 
V3 V2 V2 V3 



-1 



respectively. The flavor-charge 4x4 matrices for the 
quark-photon, diquark photon and exchange diagrams 
read: 

X>,Qti, E^TrjshQ}, E^ Sj Q s t tj , 



»=i 



where Q = diag(q u , qd) is the quark charge matrix. Sim- 
ilar expressions arise for the seagulls, see App. A. 9 of 
Ref. [47] for details. As expected from isospin symmetry, 
the flavor-charge matrices for all diagrams are propor- 
tional to diag(2, 1,0, — 1), with a coefficient 2 / 3 f° r the 
diquark part of the impulse approximation and Y 3 f° r 
all other diagrams. Consequently, the full Dirac-color- 
flavor-charge current for each member of the A multiplet 
with charge e r is given by 



j _ZL[ jIMP-Q _|_ 2 jIMP-DQ _ jEX _ jSG _ jSG 

3 



and the four different A states only differ by their charge. 
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